Conformal invariance in three-dimensional rotating turbulence 
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We examine three-dimensional turbulent flows in the presence of solid-body rotation and helical 
forcing in the framework of stochastic Schramm-Lowner evolution curves (SLE). The data stems 
from a run on a grid of 1536"^ points, with Reynolds and Rossby numbers of respectively 5100 and 
0.06. We average the parallel component of the vorticity in the direction parallel to that of rotation, 
and examine the resulting {i^z)j, field for scaling properties of its zero-value contours. We find for the 
first time for three-dimensional fiuid turbulence evidence of nodal curves being conformal invariant, 
belonging to a SLE class with associated Brownian diffusivity k = 3.6 ±0.1. SLE behavior is related 
to the self-similarity of the direct cascade of energy to small scales in this fiow, and to the partial 
bi-dimensionalization of the flow because of rotation. We recover the value of k with a heuristic 
argument and show that this value is consistent with several non-trivial SLE predictions. 
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Self-similarity in physics is a common phenomenon, 
with identical properties of a system when considered at 
different scales. Rugged coast lines, fractals, traffic in 
computer networks, growth processes, geometrical prop- 
erties of interfaces, phase transitions in critical phenom- 
ena such as in the Ising model for spontaneous magneti- 
zation, classical and quantum field theory, often display 
power-law scaling of some variable and such scaling ex- 
ponents have been the object of intense investigations 
resulting in the finding of broad classes of universality. 

A property stronger than scale invariance is conformal 
invariance, under transformations that preserve angles 
with rescaling that depends on position; it is difficult to 
test, since it implies the need to investigate the scaling 
of multi-point high-order correlation functions. However, 
recent developments by Schramm in particular (see e.g. 
[T] and references therein) allow in some cases for a sta- 
tistical characterisation of conformal invariance. Such 
scaling laws can be related to Brownian motion (which is 
scale invariant, and conformal in two dimensions) in what 
is now named Schramm-Lowner evolution (SLE), with 
as sole parameter the diffusivity k associated with this 
Brownian motion. In this approach, the driving of the 
Lowner equation (Eq. [2] below) is stochastic, with a con- 
formal map allowing to go from static (fixed time) two- 
dimensional (2D) paths in the complex plane C to "dy- 
namic" one-dimensional (ID) motions. In other words, it 
allows one to describe paths in C by a succession (convo- 
lution) of conformal maps obeying a differential equation. 
Schramm's theorem (see, e.g., |T]) states that if and only 
if the driving is Brownian is the measure of the 2D paths 
conformally invariant. 

Two-dimensional turbulence differs in many ways from 
the three-dimensional (3D) case because of the presence 
of an extra invariant in the absence of viscosity, the en- 
strophy S = (|V x up/2), leading to an inverse cascade 
of energy E = (|up/2) to large scales [5], with u the 



velocity field. It was shown in [5] that this inverse cas- 
cade, which is known to lack intermittency and is self- 
similar, can be viewed in the framework of conformal in- 
variance when examining zero-vorticity lines; it belongs 
to the universality class with k = 8/3 (the enstrophy 
cascade to small scales, however, is not SLE [3]). These 
results stem from direct numerical simulations (DNS) on 
grids of up to 16, 384^ points, with forcing at intermedi- 
ate wavenumber, kp/kmin ~ 100, with kmin — Stt/Lq, 
Lq being the size of the vessel. 

In the case of 3D Navier-Stokes (NS) incompressible 
flows at high Reynolds numbers, the cascade of energy to 
small scales is not self-similar, because of the presence of 
strong vorticity gradients. Only one time scale is present, 
the eddy turn-over time tnl ~ ^/w^, with Ui the veloc- 
ity at scale £, and dimensional analysis gives an energy 
spectrum E{k) cx fc^^/'^ that is quite close to observed 
spectra in the atmosphere or in laboratory experiments. 
However, when introducing solid body rotation fl with 
inertial time tq ~ 1/fi, E{k) steepens and its spectral 
index can be recovered by taking into account the weak- 
ening of nonlinear interactions due to the inertial waves 
[1]. In this case, self-similarity and Gaussianity in the 
3D direct energy cascade was found recently both in the 
laboratory [3 E] and in DNS [3 [5] , more clearly so in the 
presence of helicity, i.e., velocity- vorticity correlations [S]. 

Since rotating flows tend to become quasi-2D (but not 
strictly 2D, as our results will confirm) when strong ro- 
tation is imposed, the question thus arises as to whether 
SLE can be identified in such fiows. To this end, we ex- 
amine the large data set produced in a run of rotating 
helical turbulence on a grid of 1536'^ points, with = 2tt 
and forcing at kp=7] an inverse cascade of energy to large 
scales (with constant negative flux) is observed, but with 
too little extent in wavenumber to allow for a SLE analy- 
sis similar to that performed in [3] for the 2D NS inverse 
cascade. We concentrate instead on the direct energy 



2 




FIG. 1: [Color online) Perspective volume rendering of vor- 
ticity intensity in a snapshot of the flow. The slice in the 
middle of the box shows a 2D cut of ui^, while the slice at the 
bottom shows the same field averaged vertically. The inset 
is (i^z)^, with super-imposed nodal paths, the traversing ones 
(top to bottom, red online) being discarded from the analysis. 



cascade to small scales (with constant and positive flux, 
see Fig. 9 in 8J), and analyzed for its classical statis- 
tical properties and structures in [S]. A pseudo-spectral 
code with periodic boundary conditions was used, with 
at the onset of the inverse cascade a Reynolds number 
Re — Uo'in /[vkp] ~ 5100 (with v the viscosity), and the 
Rossby number Ro — U^kp /[2-kVI] k, 0.06; Uq ~ \ \& the 
r.m.s. velocity. We integrated the 3D NS equations in 
the rotating frame for an incompressible flow (V • u = 0); 
with w = V X u the vorticity, they read: 



du 

~dt 



LO X u + 20 X u : 



(1) 



V is the total pressure modified by the centrifugal term, 
and F is a helical Arn'old-Beltrami-Childress forcing [Tj 
|S]. The rotation is imposed in the vertical (z) direction, 
with ft = 9z. The code is fully parallelized, uses the 2/3 
de-aliasing rule, and the temporal scheme is a second- 
order Runge-Kutta. Note that in 3D, besides energy, 
total helicity H = {u ■ u)) is also an ideal invariant (TJ. 

The procedure: Considering the symmetries of Eq. 
we construct a 2D field by averaging in the vertical di- 
rection the parallel component of vorticity, which we de- 
note hereafter (wz)^; we also compute a transverse aver- 
age (wz)y to compare with. Starting from an arbitrary 
line, say the x-axis, we explore iso-contours of zero field 
as trajectories in the 2D plane that keep positive field 
to their right. The direction along the trajectories is 



parametrized by a "driving time" t. The path is stopped 
whenever it returns to the initial axis. The end-point 
is then sent to infinity through a holomorphic (Mobius) 
transformation as in [S], with a cut-off A chosen to be 
such that the tip of the curve is within a small arbitrary 
distance of the chosen axis; results are insensitive to the 
choice of A in a range of 1 to 10 pixels, and agree as 
well with a procedure in which the Mobius conformal 
mapping is not applied. Note that, because of periodic- 
ity, the procedure is not affected by the boundaries, and 
that all trajectories are renormalized to Tmax = 1- 

We have shown numerically for this flow the existence 
of scale invariance for the direct energy cascade and the 
Gaussianity of the velocity in [5] (see Figs. 7 and 8), also 
examining anisotropy at different times using a SO {2) x M 
decomposition (see Figs. 2 and 3 in [5] for the actual 
scaling ranges). We now probe the conformal invariance 
of these 2D curves viewed as paths in the upper com- 
plex plane; the paths are encoded in a "driving function" 
£_{t) obtained through the chordal Lowner equation be- 
low, with griC) (C € C) a conformal map (see, e.g., |12l): 



9r5r(C) 



5r(C)"^(T) 



(2) 



^(r) is the unknown ID real continuous stochastic driv- 
ing function for the path. In order to estimate ^(r) nu- 
merically, we use the zipper algorithm (ZA) with vertical 
slits rni. Then ga.SriO = a+ ^/[(C - a)^ + 4(5t] confor- 
mally maps the upper plane minus the vertical slit in C, 
[(a, 0); (a, 2^/T)], into the upper plane: ZA gradually zips 
the whole path onto the x-axis using the composition of 
functions gaMiO for different 5t. We thus transform the 
erratic nodal line in the plane (inset in Fig. [l] described 
below) into an unknown motion along the real axis, £,{t). 

To test for conformal invariance, we therefore must 
ask: Is ^(r) a Gaussian process? Does it correspond to 
a Brownian motion? And if so, what is its diffusivity? 
To answer the first question, one can use the classical 
Kolmogorov-Smirnov (KS) test, and check (i) whether 
its pks value is above a given threshold for a wide range 
of driving times r, and (ii) whether the steps in this mo- 
tion are independent. When both tests are favorable, we 
then consider the scaling of the variance of ^(r). If the 
scaling is reasonably linear with t, we will conclude that 
the set of driving functions likely stems from a Brownian 
process, and hence that the vorticity isolines obtained as 
indicated above are likely to be conformally invariant. 
The linear scaling also gives us the diffusivity k which 
describes entirely the statistics of the SLE process. 

Results: We now apply the procedure to the 1536^ 
DNS data. After performing the average (either in z or 
in y), fifteen temporal snapshots are analyzed, separated 
by approximately one eddy turn-over time. The resulting 
dataset has in excess of 3.5 x 10^ points for each aver- 
aging direction. In Fig. IT] is given a snapshot of |<^(x)| 
in 3D, a 2D slice of Wz(x) (in the middle of the box), 
and the same field component when vertically averaged, 
(wz)^ (bottom slice). The flow displays features of both 
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FIG. 2: {Color online) SLE analysis for (a;z)2 (circles) and 
(wz)y (triangles). Filled symbols comprise the full data, and 
open symbols correspond to paths in the dataset filtered at 
I < 2n/10. Driving time r is in log coordinates. Top: KS 
test; inset: PDFs of driving function increments for {oJz}^ for 
Ti and T2 as marked on the r-axis; the dashed line indicates 
a Gaussian. Bottom: Diffusivity k for (ijJz)^ with error bars; 
inset: scaling of variance for {oj^)^, with dash line indicating 
linear variation with r. 



2D and 3D behavior [7l |9] ; prominent are tangles of vor- 
tex filaments loosely organized parallel to fiz, and Bel- 
trami core vortices, smooth long-lived helical columns. 
The inset shovifs (wz)^ face-on with a few examples of su- 
perimposed paths that are analyzed below, noting that 
we discard the traversing paths, only keeping returning 
paths (loops at the top, blue online), as done in [51 [TU]. 

Figure [2] summarizes the analysis, for averaged ei- 
ther parallel (full circle) or transverse (y, full triangle) to 
fiz. Since a barely-resolved inverse cascade of energy de- 
velops in the DNS between the forcing scale and the box 
size, we also performed the analysis in a dataset in which 
was filtered so as to only preserve scales smaller than 
the driving scale: for £ < 27r/10, the results are now dis- 
played with open symbols. Figure [2] (top) gives the pks 
values of the KS test with abscissa r in log scale. The 
value Pks ~ 10"^ for (wz)^ (triangles) shows that the 
transverse y-averaged field is not Gaussian, and we shall 
not analyze further such y-averaged data. On the other 
hand, pks > 0.6 for most values of r for (a;z)z (circles). 
These opposite results imply that our test can eliminate 
non-Gaussian behavior, and that, due to the anisotropy 
of the flow, only parallel z-averaging may lead to confor- 
mal invariance. To confirm the Gaussianity of the process 
with parallel averaging, we show in the inset the proba- 




FIG. 3: {Color online) Winding angle probability Pjc/t as a 
function of the angle 9, as illustrated by the sketch at the 
bottom. Inset: Mean gyration radius D of nodal lines as a 
function of the number of pixels A'^. Dashed lines: theoretical 
predictions for ks = 3.6; open and filled circles as in Fig. [2] 



bility distribution functions (PDFs) for two driving times 
Ti and T2', the dotted line is a Gaussian with zero mean 
and unit variance |13j. Note that Gaussianity also im- 
plies independence of increments (we show the evolution 
of the PDFs for different r in Fig. [2]). 

Figure|2] (bottom) gives the variance of the motion nor- 
malized by r with error bars for the datasets that are not 
discarded by the KS test. The resulting diffusivity for the 
associated Brownian motion is (^^) jr — n — 3.6 ± 0.1 
for the full data set (full circles), and k — 3.5±0.2 for the 
data in which the scales comparable to or larger than the 
forcing have been filtered out (open circles). The inset 
shows the actual scaling of variance with driving time, 
in log-log coordinates, with the dashed line for linear 
variation. Note that, within error bars, the results are 
insensitive to whether or not we filter the numerical data 
(keeping only Fourier modes k > 10, thereby making sure 
we restrict the data to the direct cascade of energy). 

Finally, we confirm the scaling we found for k by ex- 
amining some of the predictions on statistical properties 
of nodal lines of (clIz)z that can be made using the SLE 
framework (see [T|). A classical one concerns the fractal 
dimension of the nodal lines, but less trivial features pre- 
dicted by SLE include for instance the so-called "winding 
angle" or the gyration radius. The winding angle predic- 
tion states that the probability Pieft of a SLE line to leave 
a point zq = pe*^ in C to its left depends only on k and 
B following a known expression [T] . Figure [3] shows the 
results obtained from our datasets as a function of 0, as 
well as the mean gyration radius of the nodal lines as a 
function of their length in pixels in the top-right inset. In 
both cases, the SLE predictions for k = 3.6, given with 
the dashed lines, appears convincing. 

We thus conclude that our analysis identifies confor- 
mal invariance for nodal lines of the vertical component 
of the vorticity field when averaging parallel to the di- 
rection of rotation, and fails to identify such invariance 
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in its transverse average. For the parallel-averaged ver- 
tical vorticity, the associated difFusivity is k w 3.6 ±0.1. 
Moreover, SLE predictions for this value of k agree well 
with our results. It is also important to remark that our 
analysis could fail to reject the hypothesis of Gaussianity 
if data were insufficiently resolved; this is not surprising 
since it is hard to distinguish SLE behavior from some- 
thing close to SLE [T3]. In spite of these limitations, the 
data analyzed here up to the spatial resolution considered 
is found to be consistent with SLE behavior. 

Discussion: Rotating helical turbulence may be per- 
haps the first documented case presenting SLE scaling 
for three-dimensional flows undergoing a direct cascade 
of energy and of helicity to small scales, when properly 
averaged in the direction of rotation. Conformal invari- 
ance is a strong local property and allows determination 
of a series of scaling laws, as exemplified in |3l [10] for 
2D NS and other related 2D cases such as surface quasi- 
geostrophic (SQG) flows, and as found here as well. SLE 
obtains convincingly for the vertical component of the 
vorticity averaged along the direction of rotation, with 
K « 3.6, close (but not identical) to the value identified 
in [TU] for SQG flows for an inverse cascade. Note that 
anisotropy of this 3D rotating flow must play an essential 
role, since the direct cascade of enstrophy in strictly 2D 
NS is not SLE as shown in [SJ. 

The fractal dimension D^ < 2 of SLE curves can be re- 
lated to K [3], as well as to the cancellation exponent nc 
which measures how fast a field changes sign [T5] . With d 
the dimension of space, we have Dp = 1 + k/8 = d — 2Kc- 
It is straightforward to relate the diffusivity of the SLE 
process k and the exponent e of the energy spectrum, 
E{k) ~ k~'^, under the assumption of self-similarity, — 
asp] Cp ^-re the exponents of the pth-order longitudinal 
structure functions of the velocity field, {SuLirY) ~ r'^'' : 
where 5ul is the variation of the velocity projected along 
the direction of the spatial increment r. We use that 
Kc = Ci = o,Si C2 ~ s ~ 1 for 1 < e < 3, and that dimen- 



sional analysis for a given dynamics gives 05. Then, 

'l=d-e = d- l-2as. (3) 
8 

Hence, the value of k is quite sensitive to e or as |16j . For 
2D NS, 05=1/3 and k=8/3, as found in [3] (with dual 
value K* = 6). For rotating helical turbulence, as=3/4, 
using a phenomenological model based on three assump- 
tions [3 [H] : wave- modulated energy spectrum; domina- 
tion of the helicity cascade to small scales; and maximal 
helicity. The first hypothesis allows to write that the 
transfer of energy to small scale is slowed down in the 
proportion tnl/'T'ij! the second one stems from the fact 
that the energy undergoing an inverse cascade to large 
scale, little energy is left to feed the small scales, whereas 
helicity only possesses a small-scale cascade and thus 
is the determining factor in this direct cascade. These 
two concepts lead to e + h = A, with helicity spectrum 
H{k) ~ k~'^. The third assumption gives h = e — 1 and 
thus = 3p/4, a value reported experimentally as well 
P[TS]. From Eq. we then obtain k = 4, close to the 
value we find given the statistics. 

The connection between SLE and statistical proper- 
ties of turbulence allows one to look at such flows with 
a new eye, and to build bridges between fluid dynam- 
ics and other research areas in mathematics, condensed 
matter, percolation, and quantum field theory. Other 
three-dimensional flows may be studied with the same 
tools when the flow is self-similar and symmetries allow 
for a reduction of dimensionality. As an example, we 
leave for future work an investigation of SLE properties 
in the inverse cascade of rotating turbulence. 
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